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We present results of a high statistics study of /^r, fx: fo-, /d^, and fy^ in the 
quenched approximation using Wilson fermions at /3 = 6.0 on 32"^ x 64 lattices. We find that 
the various sources of systematic errors (due to setting the quark masses, renormalization 
constant, and lattice scale) are now larger than the statistical errors. Our best estimates, 
without extrapolation to the continuum limit, are /^r = 134(4) MeV, fx = 159(3) MeV, 
/d = 229(7) MeV, fo, = 260(4) MeV, and fv^irup) = 0.33(1), where only statistical 
errors have been shown. We discuss the extrapolation to the continuum limit by combining 
our data with those from other collaborations. 
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1. Introduction 



Phenomenologically, /d, fo^, fs and fs^ are essential ingredients needed to deter- 
mine the less well know elements of the Cabibbo-Kobayashi-Maskawa mixing matrix. As 
these heavy-light decay constants are at best very poorly measured, there has been a large 
effort by many lattice groups to estimate them from numerical simulations. Decay con- 
stants are amongst the most precise quantities that one can calculate on the lattice and a 
recent review has been presented by Allton at LATTICE95 0. In this paper we present 
results for f^, /k, fx/fn^fD, fo/fn, fos, fnj fn, and vector decay constant fy^ from 
simulations done on 170 32"^ x 64 quenched lattices at /3 = 6.0 using Wilson fermions. We 
emphasize that extrapolations to the continuum limit, incorporating the results from other 
collaborations, are not very reliable as the combined data do not show an unambiguous 
pattern of 0(a) corrections. 

Preliminary results from a subset of 100 lattices were presented at the LATTICE94 
meeting f^. The raw lattice results have not changed significantly since then, however we 
now present a more detailed analysis of the systematic errors. We estimate the uncertainty 
in the results due to extrapolation of the lattice data to physical values of the quark masses, 
the renormalization constants for the lattice currents, and the extraction of the lattice scale. 
We find that these various systematic errors are now much larger than the statistical errors. 
Finite size errors, if present, are smaller than the statistical errors. Our best estimates are 
now given in the scheme called TADl to evaluate the renormalization constants for the 
axial and vector currents. 

The details of the lattices and the calculation of the spectrum are given in a com- 
panion paper . In section |] we briefly summarize the lattice parameters, and in section 
^ we describe the lattice methodology and the consistency checks made to extract the 
decay constants using estimates from different types of fits and interpolating operators. 
The choice of renormalization constants for the axial and vector currents, Za and Zy, 
is discussed in section ^, the lattice scale in section ^, and the quark masses in section 
p. In section ^ we compare the data with predictions of quenched chiral perturbation 
theory. The extrapolation of the data to physical values of quark masses is discussed in 
section Rl and our best estimates at /3 = 6.0 are summarized in section p|. In section p!0| we 



compare our data with those from other collaborations ( GFll [Q, JLQCD [Q, and APE 
1^. The MILC data presented in are preliminary and therefore not included in this 
analysis.) and extrapolate the combined data to the continuum limit. Finally, we present 
our conclusions in section |ll]. 



2. LATTICE PARAMETERS 

The details of the 170 32'^ x 64 gauge lattices used in this analysis are given in 0. We 
refer the interested reader to it for further details of the signal in the 2-point correlation 
functions and on the extraction of the spectrum. In this paper we concentrate on the 
analysis of systematic errors in decay constants associated with fixing the quark masses 
^ = {i^u + "id)/2, rris and rric, the renormalization constants Za and Zy, and the 
extrapolation to physical masses and the continuum limit. 
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To calculate decay constants we used the Wuppertal source quark propagators at five 
values of quark mass given by a? = 0.135 (C), 0.153 (5), 0.155 (t/i), 0.1558 (t/2), and 0.1563 
(C/3). These quarks correspond to pseudoscalar mesons of mass 2835, 983, 690, 545 and 431 
MeV respectively where we have used 1/a = 2.33 GeV for the lattice scale. We construct 
two types of correlation functions, smeared-local (Tsl) and smeared-smeared (Fg^) which 
are combined in different ways to extract the decay constants as discussed below. The three 
light quarks allow us to extrapolate the data to the physical isospin symmetric light quark 
mass m, while the C and 5" k values are selected to be close to the physical charm and 
strange quark masses. The physical value of strange quark lies between 5" and Ui and we 
use these two points to interpolate to it. In most cases we find that extrapolation to fn can 
be done using the six combinations of light quarks C/if/i, f/if/25 UiU^, U2U2, U2US, UsUs. 
For brevity we will denote this combination by {UiUj} and the three degenerate cases by 



3. Lattice Method for calculating fps and fy 

The lattice definition of the pseudo-scalar decay constant fps, using the convention 
that the experimental value is /^r = 131 MeV, is [§] 

_ ZA{0\Ai-^^'\n{p)) 

where Za is the axial current renormalization constant connecting the lattice scheme to 
continuum MS. In order to extract /tt we study, in addition to the 2-point correlation 
functions F, two kinds of ratios of correlators: 



TsLit) t^oo {0\A,'''''^'\n) 



Ri{t) 
R2{t) 



Tssit) 2M^ 



In the case of Ri we have to extract (01^4*^™'^^'^'^'^ |7r) separately from the F55 correlator. 
For each of the two ratios, Ri and -R25 the smeared source J used to create the pion can 
be either tt or A4. This gives four ways of calculating f^^, which we label as (using ratio 
Ri with J = tt), (using ratio Ri with J = A4), (using ratio R2 with J = tt), and 
(using ratio R2 with J = A4). The fifth way, consists of combining the mass and 
amplitude of the 2-point correlation functions {A4P)ls and {PP)sSj and the sixth way, 
fl, uses (^4^4)1,5 and (^4^4)55. 

The lattice results for mesons at p = for the different combinations of quarks are 
given in Table |l| using the renormalization scheme Ztadi defined in Table |^. All Errors 
are estimated by a single elimination Jackknife procedure. The results from the six ways 
of combining the two correlators are mutually consistent. Since the different methods use 
the same correlators, the data are highly correlated; however, consistent results do indicate 
that fits have been made to the lowest state in each of these correlators and reassure us of 
the statistical quality of the data. 
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Table 1. The data, in lattice units, for the pseudo-scalar decay constant fps for the six different 
ways of combining the SL and SS correlators described in the text. The renornial- 
ization scheme used to generate this data is TADl as described in Table and the 
meson mass used in the analysis is taken to be the pole mass. 



The results for fj^ using correlators at non-zero momentum are given in Table ^. The 
data show that in almost all cases the results are consistent within 2a. The most noticeable 
differences are in the p = (2, 0, 0) values for lighter quarks. The signal in these channels is 
not very good and it is likely that in these cases there exists contamination from excited 
states over the range of time-slices to which fits have been made. We regard the overall 
consistency of the data as another successful check of the lattice methodology. Henceforth 
we shall restrict the analysis to p = (0, 0, 0) case as it has the best signal. 

The dimensionless vector decay constants are defined as 

Zv{0\V^'°''''\V) = (3.3) 

Jv 

where is the vector current and \V) is the lowest 1~ state with mass My. The experi- 
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Table 2. The data, in lattice units, for the pseudo-scalar decay constant fps, averaged over 
the six different ways of combining the SL and SS correlators, measured at different 
momenta. The renormalization scheme is TADl as described in Table ^ and the 
meson mass used in the analysis is taken to be the pole mass. 

mental quantities are related to fy^ by 

-^fv\M,) = 0.199(5), 

lfy\M^) = -0.078(1), (3.4) 

p-\Mj/^) = 0.087(3), 

where the values are calculated using the rate r{V — > e+e") given in the PDB94 Q. We 
extract the relevant matrix element in the same two ways as described in Eq. ( |3.2|) for fps- 
To study discretization errors we study three lattice transcriptions of the vector current 



J p 

f(j} 

•I j/ip 
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(local, extended, and conserved). 



V^{x) = i^{x)-ff,Uf,{x)ij{x + ij)+i^{x + il)-f^Ul{x)tP{x), (3.5) 
= ij{x){'y^,-r)Uf,{x)ij{x + fi) +i){x + fl){'yf, + r) Ul{x)ij{x). 

where for degenerate quarks the last form is the conserved current. In Tables ^ and 
^ we show the lattice data for the 15 mass combinations as a function of the different 
methods/currents, and versus the renormalization schemes for the local current. Overall, 
the data show that the two methods in Eq. ( |3.2| ) give consistent results for all three currents. 
The results from the local and extended vector currents also agree, while those from the 
conserved current are ~ 10% smaller. These points will be discussed in more detail later. 

In order to extract results that can be compared with experiments we analyze the 
data in terms of the five sources of systematic errors discussed below. 



4. The renormalization constant Za and Zy 

Reliable calculations of decay constants depend on our ability to calculate the renor- 
malization constants, Za and Zy, linking the lattice and continuum regularization schemes. 
In our analysis we use 1-loop matching with the tadpole subtraction scheme of Lepage- 
Mackenzie. An outline of the scheme, which includes picking a good definition of the lattice 
as and the scale q* at which to evaluate it is as follows. Lepage and Mackenzie show that 
(to be defined below) is a better expansion parameter than the bare lattice coupling. To 
pick the value of g* we need to know the "mean" momentum fiow relevant to a given matrix 
element. Again it has been pointed out by Lepage and Mackenzie that q* , estimated by 
calculating the mean momentum in the loop integrals, is dominated by tadpole diagrams 
which are lattice artifacts. If the tadpoles are not removed then this scale is typically n/a. 
They have proposed a meanfield improved version of the lattice theory which removes the 
contribution of tadpoles. The effect of this is three-fold. One, it typically changes q* to 
1/a, i.e. the matching scale becomes more infrared if the tadpole diagram is removed; sec- 
ond the renormalization of the quark field changes from \/2k — > a/8kca/1 — 3k/4kc', and 
finally the perturbative expression for 8kc is combined with the coefficient of in the one 
loop matching relations to remove the tadpole contribution. 

To get Q!j\^g(g*) we use the following Lepage-Mackenzie scheme. The coupling cxy is 
defined at scale q = 3.41/a to be 

ay(^^) (1- (1.19 + 0.017nf)ay) = -—ln(-TrPlaq), (4.1) 
a 47r 3 

which is related to aj^g at scale q = 3.41/a by 

^— = — + 0.822. (4.2) 
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Table 3. Lattice data for the vector decay constant fy^ for the two different ways of combining 
the SL and SS correlators, and for the three different lattice vector currents described 
in the text. The renormahzation scheme in all cases is TADl as described in Table ^ 
and the meson mass used in the analysis is taken to be the pole mass. 



We then run o;j^ from q to q* by integrating the 2-loop /3-function. To translate the 
results from q* to any other point one uses the standard continuum running. 

At the lowest order there are two equally good tadpole factors, Uq = plaquette^^^ or 
8kc- To the accuracy of the meanfield improvement one expects SkcUq = 1. Deviations 
from this relation (^ 10% for the Wilson action at /3 = 6.0) are a measure of possible 
residual errors. Writing the tadpole factor as 1 — Xaj^jg{q*), we define a given Z factor 
to be 

Z = 1 + aj^iql log(Q*a) + [C - X)^ (4.3) 
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Table 4. Lattice data for the vector decay constant fy^ as a function of the different renornial- 
ization schemes given in Table |. The results are for the local current, and the meson 
mass used in the analysis is taken to be the pole mass. 



where C is the difference between the finite part of the continuum MS and lattice 1-loop 
result. Thus, Za for the local operator in the tadpole improved schemes is 

^/Z^Z^ = Vl - 3a^i/4k,V1 - 3a^2/4kc (1 - a^(?*)(1.68 - X)). (4.4) 

In order to examine the dependence of the decay constants on Z and the renormaliza- 
tion of the quark field we present our results for seven different commonly used schemes 
described in Table |^. The schemes ZxADa, Ztadi, Ztad2j Ztad-k^ and Zxgf 11 are all 
self-consistent to 0{as)- The scheme ZtadUo ^^'^ have replaced 8kc by Uq in 

only one part. We shall quote, as our best estimates, results obtained in the Ztadi scheme 
and use the difference between it and Ztad-k as an estimate of the systematic error due 
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to tuning q* . Finally, an estimate of the residual perturbative errors is taken to be the 



difference between Ztadi and ZtadUq-: and is given in column labeled Za in Table ^ 



This, we believe, is an over-estimate of the error we make by using the 1-loop coefficient 
of a,,. 
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Tv/a 




0.202 


0.190 


0.151 


0.133 


0.190 
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0.133 



Table 5. The different renormalization sdiemes used in tlie analysis. Tlie two possible tadpole 
factors are Uq = plaq^^^ = 0.878 and 1/8kc = 0.795. The 1-loop perturbative expan- 
sions for these are Uq = 1 — 1.0492q! and 8kc = 1 + 1.364q: respectively. The sixth 
scheme ZTgfu is the one used by the GFll collaboration with a slightly different 
definition of ctj^^. 



The renormalization of the local vector current, Zy proceeds in the same way as Za- 
In case of both the extended and conserved currents, there is no tadpole contribution 
in C as it cancels between the wave-function renormalization and the vertex correction. 
Consequently, we use the non-perturbative value for 8kc. The complete renormalization 
factors in the tadpole improved schemes for relating the lattice results to the continuum 
are 

^ZlZlZ^ = Vl - 3kJAk,^1 - 3k2/Ak, (l - (2.182 - X)), 

^Jz^ZlZ^ = 8KeVl - 3ki/4kcV^ - 3k2/4k, (l - 1.038aj^{q*)), (4.5) 

^Z^Z^Z^ = 8k,^1 - 3ki/4k,^1 - 3k2/4k,. 

We find that the results with the local current lie in between those from the extended and 
conserved currents, and have the best statistical signal. We therefore quote results from 
the local current as our best estimate, and use the difference between them as an estimate 
of the systematic error. 



5. The lattice scale a 

To convert lattice results to physical units we use the lattice scale extracted by setting 
Mp to its physical value. This gives 1/a = 2.330(41) GeV . The variation of 1/a between 
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the Jackknife samples is folded into our error analysis, however different ways of setting 
the scale are not. For example, using to set the scale gives 1/a = 2.018(37) GeV 0, 
while NRQCD simulations of the charmonium and T spectrum give 1/a = 2.4(1) GeV 
||10|| . As we show later, the scale determined from f-,^ is 2265(57) MeV. Thus, estimates 
based on mesonic quantities like Mp, heavy- heavy spectrum, and /tt all give consistent 
results. We take l/a(Mp) = 2.330(41) GeV and use the spread, ~ 70 MeV ~ 3%, as our 
best guess of the size of scaling violations relevant to the analysis of the decay constants. 
To reduce this error requires using an improved gauge and fermion action, which is beyond 
the scope of this work. 



6. Setting the quark masses 

In order to extrapolate the lattice data to physical values of the quark mass we have 
to fix m, rUs and nic. The chiral limit is determined by linearly extrapolating the data for 
to zero using the six cases {UiUj}. Our best estimate is 

= 0.157131(9), (6.1) 

which is used in the calculation of Z^. 

To fix the value of ki corresponding to m we extrapolate the ratio M^/M^ to its 
physical value 0.03182. The result is 

Ki = 0.157046(9) (6.2) 

Thus, our data is able to resolve between the chiral limit and m. In we had shown that 
a non-perturbative estimate of quark mass m„p, calculated using the Ward identity, and 
(1/2k — 1/2kc) are linearly related for light quarks, so either definition of the quark mass 
can be used for the extrapolation. We have chosen to use runp in this paper. 

The determination of the strange quark mass has significant systematic errors as 
shown below. We determine Kg in three ways as described in [Q]. We extrapolate Mj./M^, 
MK*/Mp, and M^j^/Mp to m and then interpolate in the strange quark to match their 
physical value. In Table ^ we give Kg, the non-perturbative estimate rrinpa = nisa, and 
nis = .^m(l/2K — 1/2kc) evaluated at 2 GeV in the MS scheme using the TADl matching 
between lattice and continuum. The data show a ~ 20% difference between various esti- 
mates of nis which cannot be explained away as due to statistical errors. Using M'^/M^ 
to fix nis implies that nis = '25rn as we use the lowest order chiral expansion to fit Mpg 
data. On the other hand M^/Mp gives nis/m ^ 30. This estimate is not constrained 
by the chiral expansion, and is in surprisingly good agreement with the next-to-leading 



chiral result |jTl|]. In this paper we shall quote results for both ms(Mx) and ms(M(^), 
and take the values with ms(M0) as our best estimates. The difference in results between 
using msiMx) and ms(M(^) will be taken as an estimate of the systematic error due to the 
uncertainty in setting nis- 
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K 


Til (1 


m (2 GeV) MeV 


Ml/Ml 


0.15503(7) 


0.0372(14) 


129(2) 


Mk*/Mp 


0.15479(19) 


0.0421(36) 


145(9) 


M^/Mp 


0.15464(17) 


0.0445(32) 


154(8) 



Table 6. Estimates of iris using different combinations of hadron masses. We give the k val- 
ues, the quark mass determined by the Ward identity, and = Zmi^/^K — 1/2kc) 
evaluated at 2 GeV in the MS scheme, and using the TADl tadpole subtraction 
procedure. 



To determine the value of k corresponding to mc we match Mj), Mo* and Mo^ as 
these are obtained from the same 2-point correlation functions as used to determine the 
decay constants. Unfortunately, as shown in P] , the estimate of charmonium and D meson 
masses measured from the rate of exponential fall-off of the 2-point function (pole mass or 
Ml) and those from the kinetic mass defined as M2 = {d"^ E / dp'^\p=o)~^ are significantly 
different. We find that the data for the heavy- heavy and heavy-light qq combinations are 
consistent with the nearest-neighbor symmetric-difference relativistic dispersion relation 
sinh^(E/2) - sin2(p/2) = sinh^(M/2), in which case M2, as defined above, is given by 
sinhM. The results for Mi and M2 for the D states are given in Table ^for k — 0.135 (we 
have simulated only one heavy quark mass) . The data show that the experimental results 
lie between Mi and M2 for each of the three states, and the difference between Mi and M2 
is large and statistically significant. The size of this systematic error and the uncertainty 
in setting the scale 1/a, makes it difficult to fix Hcharm- We simply assume that k = 0.135 
corresponds to and quote final results using M2. As an estimate of systematic errors 
associated with not tuning mc we take the difference in results between using Mi and M2 
since we do not have access to the rate of variation of decay constants in the vicinity of 



7. Quenched approximation 

In the last couple of years it has been pointed out by Sharpe and collaborators [0 
and by Bernard and Golterman [jl^ that there exist extra chiral logs due to the ry', which 



is also a Goldstone boson in the quenched approximation. These make the chiral limit of 



quenched quantities sick. To analyze the effects of quenching Bernard and Golterman [p!3 
have constructed the ratio 

R= (7.1) 



fivf: 



22' 



applicable in a 4-flavor theory where mi = my and m2 = m2'. The advantage of this ratio 
in comparing full and quenched theories is that it is free of ambiguities due to the cutoff 
A in loop integrals, and O(p^) terms in the chiral Lagrangian. The chiral expression for R 
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Mt^ (m (MtA^ 
^^J-Ds\''''s\-'-^J-K ) ) 






1 QfiQ 




1914(26) 


2137(27) 


1969 


MD*{ms{MK)) 


1961(31) 


2201(34) 


2110? 


MDt{ms{M^)) 


1978(27) 


2224(29) 


2110? 



Table 7. A comparison of lattice estimates of D meson masses with the experimental data. We 
show results for Mi and M2 and for the two different ways of setting ms described in 
the text. 



in the quenched theory is 



= 1 + 5 



m 



12 



log 



11' 



m: 



22' 



+ 0((mi -m2)2), 



(7.2) 



where 5 
QCD is 



ml/24n parameterizes the effects of the r]'. The analogous expression in full 



1 + 



,2 T„"^ll' 



mil, Ln 



12 



2 T „ "^22' 



+ m22' Lii 



m 



12 



+ 0((mi - 7712) 



(7.3) 



The leading analytic corrections in both cases are 0{{mi—m2)^) fl^j s-nd were not included 
in the analysis presented at LATTICE 94 |T^ . The data, shown in Figs. |1| and indicates 
the need for including them in the fits. {X quenched is the coefficient of S in Eq. ([7.2|), 
and Xfuii is the complete chiral logarithm term in Eq. (|7.3|) .) The fit to the quenched 
expression, Fig. |1|, gives 6 = 0.14(4). The fit to full QCD expression has smaller if we 
leave the intercept as a free parameter. In that case the fit gives 1.69(45) and not unity as 
required by Eq. ( |7.3| ). Thus, the effect of chiral logs is small, barely discernible from the 
statistical errors, and partly due to normal higher order terms in the chiral expansion. We 
shall therefore neglect the effects of quenched chiral logs in this study, and only discuss 
deviations of fps from a behavior linear in at the appropriate places. 

The second consequence of using the quenched approximation is that the coefficients 
in the chiral expansion are different in the quenched and full theories. This difference can 
be evaluated by comparing quenched and full QCD data, which is beyond the scope of this 
work. Thus, we cannot provide any realistic estimates of errors due to using the quenched 
approximation. 
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Fig. 1: Bernard-Golterman ratio R versus (mi — 7712)^. X quenched is the coefE- 
cient of d defined in Eq. (\7.2i). The intercept gives d — 0.14(4). 



8. Extrapolation in quark masses 

In Fig. 1^ we show the pseudoscalar data for {UiUj} and {SS,SUi} combina- 
tions along with two different linear fits, one to the six {UiUj} data points (/psa = 
0.0572(14) + 0.51(2)ma) and the other to the four SS and {SUi} points {fpsa = 
0.0568(14) + 0.48(l)ma). Here m is the average mass of the quark and anti-quark. The 
data show that even though the slopes for the two fits are different, the values after ex- 
trapolation are virtually indistinguishable. The size of the break between the {SS, SUi} 
and {UiUj} cases at mg is right at the la level, and no such break is visible between the 
UiUi and the U2U2 cases. We thus extrapolate to /^r using {UiUj} points and assume that 
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Fig. 2: Bernard-Golterman ratio R versus (mi — 1712)'^ ■ Xf^u is the chiral cor- 



rection defined in Eq. (7^). The hnear fit gives an intercept of 1.69(45) 
instead of unity as indicated by Eq. (ITT 



the overall jackknife error adequately includes the uncertainty due to extrapolation. 

In Fig. ^ we show the extrapolation for heavy-light mesons for three cases of "heavy" 
(C, 5", Ui) quarks. The linear fits in the light quark mass, 

fpsa = 0.103(3) + 0.33(5)m^pa (CU,) , 

fpsa = 0.074(1) + 0.26(l)m„pa {SU^) , (8.1) 
fpsa = 0.067(1) + 0.25(l)m„pa {UiU,) , 
fit the data extremely well in each of the three cases. Deviations from linearity are apparent 
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0.02 0.04 0.06 0.08 

Fig. 3: Plot of data for f^^a versus rrinpa. The linear fit, shown as a solid line, 
is to the six {UiUj} points. The error estimate on the fit is shown by 
the dotted lines. The dash-dot line is a linear fit to the four SS and 
{SUi} points. The vertical line at rrinpa ~ represents ffi and the band 
at 0.04 denotes the range of m^. 



if the "light" quark mass is taken to be S as shown by the fourth point at rrinpa = 0.076. 
These can be taken into account by including corrections, i.e. chiral logs and/or a quadratic 
term. A fit including a quadratic term fits all four points exceeding well, however the 
extrapolated value changes by < 0.2a in all three cases. Also, the change in curvature 
between UiUi and CUi is within the error estimates. Considering that the form of the 
correction term is not unique, and that the linear and quadratic fits give essentially the 
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same result, we consider it sufficient to use a linear fit to the three Ui points to extrapolate 
the heavy-light decay constants to m. 

The difference in slope between fits to {UiUj} and {SUi} points does effect the value 
of fx- We, therefore, calculate it in two ways; the central value is taken by extrapolating 
the {SUi} and {UiUi} data in the light quark to m and then interpolating in the "heavy" 
to rus- In the second way we use the slope determined from {UiUj} points and extrapolate 
to m + nis. The two give consistent results and we use the difference as an estimate of the 
systematic error. 
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Fig. 4: Extrapolation of heavy-light pseudoscalar decay constants for three 
cases of "heavy", C, S, Ui quarks. The linear fits are to to the three 
"light" Ui quarks, and the fourth point (light quark is S) is included to 
show the breakdown of the linear approximation. 



15 



The analogous plots for fy^ are shown in Figs. ^ and ^. To extract we make 
linear fits to the six {UiUj} and the three {UiUi} points. As shown in Fig. these two fits 
are almost identical (/\/a = 0.328(10) + 0.33{23)mnpa) and neither of them fits the data 
very well. The {SUi} points show a very significant break from the {UiUj} points, so to 
extract fx*, fn* we use the fits shown in Fig. p. As in the case of /ps, a linear fit to the 
three cases {CUi, SUi, UiUi) works well. The fit parameters are 

/v^a = 0.163(6) + 0.31(10) m^pa (CUi) , 

fva = 0.300(9) + 0.026(14)m„pa {SU^) , (8.2) 

/i/a = 0.322(7) - 0.18(10) m^^a (UiU,) , 

Note that the slope changes sign between the SUi and UiUi cases. Since the points at 
TfinpCL = 0.076 {S) show deviations from the linear fits, we do not include this point in our 
analysis. 



9. Results at /3 = 6.0 

The results for the pseudoscalar decay constants, in lattice units, are given in Table ^ 
for each of the seven renormalization schemes. The table also shows the variation with 
respect to the two choices of and whether one uses Mi or M2 for the heavy-light meson 
mass. For our best estimates we use Ztadi^ and convert this data to MeV using l/a{Mp). 
The results are summarized in Table ^ where we again display variation with respect to 
rus and the heavy-light meson mass. 

Our final results are shown in Table |l^ along with the estimates of the various sys- 
tematic errors discussed above. Thus, at /3 = 6.0 the value of /tt come out about 3% 
larger. Using /^r data to set the lattice scale gives l/a(/7r) = 2265(57) MeV, whereas 
l/a{Mp) = 2330(41) MeV Even ignoring the various systematic errors, the two esti- 
mates differ by roughly la. 

The ratio JkI J-k = 1.186(16) is about 2a smaller than the experimental value 1.223 if 
one ignores all systematic errors. (The systematic error in fixing mg would tend to lower 
our estimate, i.e. further increasing the difference.) An under-estimate of this ratio in the 



quenched approximation is consistent with predictions of quenched CPT ||T^ [|T3 . 

The major uncertainty in the results for the heavy-light cases, and Jdsi comes from 
the uncertainty in Zjs^ and in setting the charm mass. These corrections can be significant, 
and we need to reduce the various sources of systematic errors in order to extract reliable 
continuum estimates. 

In Tables ^ and |12| we give the values for the vector decay constant /^^, extrapo- 
lated to the masses of a number of vector states even though some of them do not decay 
electromagnetically to /"*"/". These tables also give the variation with respect to setting 
ms, the heavy-light meson mass {M\ or M2), Q*, Za-, and the dependence on the lattice 
current. The criteria that the three types of currents should give consistent results justi- 
fies using the Lepage-Mackenzie procedure for also, as pointed out by Bernard in |T^ . 
Using the ^J2k, normalization for {i.e. the same normalization as {p^^ = 0) which 
is constrained by the value of the conserved charge) gives significantly smaller values for 
cases with C quarks. 
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Fig. 5: Plot of data for fy^ versus rrij^pa. The linear fit is almost identical for 
the six {UiUj} or three {UiUi} points. The SS and {SUi} points are 
also shown for comparison. 



10. Infinite volume continuum results 

In the companion paper analyzing the meson and baryon spectrum 0] we show that 
there are no noticeable differences between results obtained on 24"^ (earlier calculations) and 
our 32'^ lattices. Thus we do not apply any finite size corrections to our data. To extract 
results valid in the continuum limit we combine our data with those from the GFll (/3 = 
5.7, 5.93, 6.17) |], JLQCD (/3 = 6.1, 6.3) [|, and APE {/3 = 6.0, 6.2) § Collaborations. We 
have attempted to correct for as many systematic differences, however some like differences 
in lattice volumes, range of quark masses analyzed, and fitting techniques remain. 
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Fig. 6: Extrapolation of heavy-light vector decay constants for three cases of 
"iieavj", C, S, Ui quarks. The linear fits are to to the three ''light" Ui 
quarks, and the fourth point (light quark is S) is included to show the 
breakdown of the linear approximation. 



We first compare the data for f^^ and fx from the different collaborations as shown 
in Figs. and ^. The various calculations have similar statistics (within a factor of two) 
and the two largest physical volumes used are by GFll (24"^ at /3 = 5.7) and LANL (32"^ 
at /3 = 6.0) collaborations. To facilitate comparison we make three changes; (a) we switch 
to the convention in which /tt = 93 MeV, (b) use the ZTgfii scheme, and (c) set nis using 
Mk. a noticeable difference in the data shown in Figs. ^ and |^ is that the APE points 
at f3 = 6.0 lie about la higher than LANL's and the value of nipa is also larger. We 
believe that the difference is partly a result of extrapolation from heavier quarks (APE 
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Table 8. Summary of results for pseudoscalar decay constants in lattice units. The variation 
with TTLs (set by Mk or Mij, ) and the heavy-light meson mass (Mi or M2 ) are shown 
explicitly. The Jackknife error estimates include statistical and a part of systematic 
errors due to extrapolation in quark masses. 



collaboration use a linear fit to extrapolate data at k = 0.153,0.1540.155 to the chiral 
limit). We find that both fps and My data show negative curvature, and a linear 
extrapolation using only SS and UiUi points increases LANL estimates, accounting for 
the full difference in Mp and a part of that in /tt . The more important feature of the data, 
however is that neither plot shows a clear a dependence. Nevertheless, a linear fit to all 
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Table 9. Results for decay constants in the Ztadi scheme as a function of iris and heavy-hght 
meson masses. The data have been converted to MeV using Mp to set the scale. Only 
the jackknife error estimates are given. 
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q* 


Tuning 
a (3%) 


Za 


u 
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4 






+2 


4 


10 


fK 


159 


3 


-3 




+3 


5 


10 


fn 


229 


7 




+12 


+4 


7 


14 


fOs 


260 


4 


-5 


+15 


+4 


8 


20 


fK/ fn 


1.19 


0.02 


-0.025 











fo/ f-n- 


1.71 


0.05 




+0.09 






? 


foj fo 


1.135 


0.021 


-0.023 


+0.006 










Table 10. Our final results using TADl scheme along with estimates of statistical and various 

systematic errors as described in the text. All dimensionful numbers are given in MeV 
with the scale set by Mp. For the systematic errors due to m^, mc, q* we also give 
the sign of the effect. We cannot estimate the uncertainty due to using the quenched 
approximation. Also, we do not have useful estimates for entries marked with a ?. 
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^TADa 






^TAD-K 


7 , 

^TADSk 


7 


7 

^ Boost-K 


(Ml) [Mp 


0.324(10) 


0.328(10) 


0.340(11) 


0.346(11) 


0.305(09) 


0.331(10) 


0.345(11) 




0.o19(0d) 


0.322(0d) 


0.335(07) 


0.340(07) 


0.300(06) 


0.325(06) 


0.331(06) 


(Mi,M^) [Mk*\ 


0.o15(0d) 


0.31o(0d) 


0.330(0d) 


0.336(06) 


0.296(06) 


0.321(06) 


0.325(06) 




0.312(04) 


0.31d(04) 


0.328(05) 


0.333(05) 


0.293(04) 


0.318(04) 


0.316(04) 




O.o0o(04) 


0.311(04) 


0.323(05) 


0.328(05) 


0.289(04) 


0.314(04) 


0.309(05) 


(Ml) [Ml)*] 


0.1d2(05) 


0.1d4(05) 


0.170(0d) 


0.173(06) 


0.152(05) 


0.165(06) 


0.134(04) 


(ikZ2) [Md* 


O.lo/ (04) 


n 10^/'^/l^ 
0.139(04) 


0.144(04) 


0.14/(04) 


n 1 on / n /I ^ 
0.129(04) 


0.140(04) 


0.114(03 ) 


(Mi,Mk) [Mz)*; 


0.173(03) 


0.175(03) 


0.182(03) 


0.185(03) 


0.163(03) 


0.177(03) 


0.140(03) 


{M2,Mk) [Md*, 


0.145(02) 


0.147(02) 


0.152(03) 


0.155(03) 


0.136(02) 


0.148(02) 


0.117(02) 




0.175(03) 


0.177(03) 


0.184(03) 


0.187(03) 


0.164(03) 


0.179(03) 


0.141(03) 


(M2,M<^) [Moj: 


0.146(02) 


0.148(02) 


0.154(03) 


0.156(03) 


0.138(02) 


0.149(02) 


0.118(02) 



Table 11. Results for fy extrapolated to the masses of a number of vector states specified 
within [] as a function of the renormalization schemes, rUg (Mk or M<^ ), and meson 
mass (Ml or M2). 



data, assuming that lattice spacing errors are 0(a), gives 

20), 

(10.1) 



^ = 0.110(5) (expt. 0.120), 



^ = 0.121(4) (expt. 0.147). 

with X^/dof = 1-6 and 1.7 respectively. The change from the GFll results is marginal as 
the fit is still strongly influenced by the point at = 5.7, which may lie outside the domain 
of validity of the linear extrapolation. A linear extrapolation excluding the P = 5.7 data 
gives 

# = 0.118(10), 



^ = 0.132(8), 



(10.2) 



with x^/dof = 2.1 and 1.9 respectively. Using ms{M^) would increase fx by ~ 2%. Given 
the large difference in the extrapolated value depending on whether the data at /? = 5.7 is 
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(Ml) [Mp] 


0.328(10) 


0.335(2d) 


0.304(18) 


[Mi,Mk) [Mk*\ 


0.o22(0d) 


r\ 000/1 0\ 

0.338(13) 


0.297(10) 


{Mi,M^) [Mk*\ 


r\ 010/ r\/^ \ 
0.dl5(0D) 


0.334(12) 


r\ ono / r\n\ 

0.293(09) 


/Tl/f 7i/f \ TTI/TI 

(Mi,Mk) [M^J 


0.310(04) 


r\ 000/ r\/^ \ 

0.332(0d) 


0.291(0d) 




0.dll(04) 


0.327(06) 


0.287(05) 


(Ml) [Md*\ 


0.1d4(05) 


0.171(0d) 


0.151(05) 


(M2) [Md*, 


O.ioy(04 ) 


0.140(00 ) 


0.128(04) 


{Mi,Mk) [Md*] 


0.175(03) 


0.182(03) 


0.163(03) 


(M2,Mk) [Md*] 


0.147(02) 


0.152(03) 


0.137(03) 


(Mi,M^) [Mz,j] 


0.177(03) 


0.183(03) 


0.164(03) 


(M2,M^) [Mz,j] 


0.148(02) 


0.153(03) 


0.138(03) 



Table 12. Results for fy as a function of the different discretizations of the vector current. We 
also show the dependence on rus (Mk or M^J and meson mass (Mi or M2). 



included or not makes it clear that more data are required to make a reliable extrapolation 
to the continuum limit. 

The and f^^ data are combined with results from JLQCD and APE ||] collab- 
orations as shown in Fig. |^. The results are in TADl scheme, and for comparison we use 
itIs^Mk)- Also, from here on we switch back to the convention in which /tt = 131 MeV. 
The APE collaboration use Mi for the meson mass. For consistency we have shifted their 
data to M2 using our estimates given in Table [1^. A linear extrapolation to a = then 
gives 

fo = 186(29) MeV, 

(10.3) 

fo, = 218(15) MeV, 

with x^/dof = 2.2 and 2.0 respectively. Using ms{M(f,) would increase fn^ to 224(16) MeV. 
The quality of the fits is, however, not very satisfactory. We feel that in order to improve 
the reliability of estimates in Eq. ( |10.3|) one needs to reduce the various systematic errors 
that have not been included in the a —>■ extrapolations presented above. 

Finally, a linear fit to data is shown in Fig. ^ The extrapolated value, 0.16(2) 
with x^/dof = 1-8, is smaller than the experimental value 0.199(5), and also smaller than 
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0.14 




I I I I \ I I I I I I I I I I I 

0.2 0.4 0.6 

Fig. 7: Linear extrapolation to the continuum limit of the ratios f-j^/Mp. Our 
data is shown with the symbol octagon, squares and fancy squares are 
the points from the GFl 1 Collaboration jj^, diamonds are APE collab- 
oration points J^, and the plus symbol labels JLQCD jj^ data. The 
two GFll points at Mpa ~ 0.56 represent 16^ (squares) and 24'^ (fancy 
squares) lattices at (3 = 5.7. 



that from a fit to just the GFll data which gives 0.18(2) g. 



11. Conclusions 

We have presented a detailed analysis of the decay constants involving light-light and 
heavy-light (up to charm) quarks. We find that the various sources of systematic errors 
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0.2 



0.4 



0.6 



Fig. 8: Linear extrapolation to the continuum limit of the ratios fx/Mp. Our 
data is shown with the symbol octagon, the squares and fancy squares 
are the points from the GFl 1 Collaboration and the diamond labels 
APE §/ data. 



(due to setting the quark masses, renormalization constant, and lattice scale) are now 
larger than the statistical errors. Work is under progress to address these issues. Our best 
estimates for the pseudo-scalar decay constants and the various sources of error, without 



extrapolation to the continuum limit, are given in Table 10 



We would like to stress that including all of the present high-statistics large lattice 
data, the extrapolation to the continuum limit is, in all cases, not very reliable. For the 
Wilson action the corrections are 0(a), and one expects that a linear extrapolation should 
suffice starting at some /3. We find that in all cases the combined world data do not show 
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IQQ I I I I \ \ I I I I I I I I I I I I I I I 

0.1 0.2 0.3 0.4 

Mp a 

Fig. 9: Extrapolation to the continuum limit of fo and Jd^ (in MeV) data. 

Our data is shown with the symbol octagon, the plus points are from 
the JLQCD Collaboration and the diamonds label the APE collab- 
oration ^ data. 

an unambiguous linear behavior in a. Since different groups analyze the data in different 
ways, there is no clean way of including the systematic errors in individual points in the 
fits. We, therefore, cannot resolve whether the poor quality of the linear fits is due to the 
various systematic and statistical errors or due to the presence of higher order corrections. 
As a result, our overall conclusion is that precise data at a few more values of jS are required 
in order to extract reliable results in the a —>■ limit. 

We have made linear fits to the data with and without including the point at the 
strongest coupling, /3 = 5.7. A linear fit to combined world data gives /tt = 120(6) MeV 
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0.2 



J I I L 
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J I I L 
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0.6 



Fig. 10: Extrapolation to the continuum limit of fp^- Our data is shown with 
the symbol octagon, and the rest of the points are from the GFll Col- 
laboration /H/. 



and fx = 135(5) MeV. Excluding (3 = 5.7 point changes these estimates to /tt = 128(6) 
and fx = 146(5) MeV. Our best estimates for heavy-hght meson, /d = 186(29) MeV and 
/x)^ = 218(15) MeV in the continuum hmit, are from a hnear fit to data at /3 > 6.0. The 
above estimates are using ms(Mx). Using ms(M^) (our preferred value) would increase 
fx and fo, hy ^ 2%. 

We study three lattice transcriptions of the vector current to calculate fy^. Using 
the Lepage-Mackenzie scheme to calculate Zy for each of the three currents yields results 
that are consistent to within 10%. We extrapolate to the continuum limit by com- 
bining with results from the GFll collaboration. The result is 0.16(2) compared to the 
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experimental value of 0.199(5). 
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